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Abstract. A microscopic approach to the investigation of the behaviour of a 
symmetrical binary fluid mixture in the vicinity of the vapour-hquid critical point 
is proposed. It is shown that the problem can be reduced to the calculation of 
the partition function of a 3D Ising model in an external field. For a square-well 
symmetrical binary mixture we calculate the parameters of the critical point as 
functions of the microscopic parameter r measuring the relative strength of interactions 
between the particles of dissimilar and similar species. The calculations are performed 
at intermediate (A = 1.5) and moderately long (A = 2.0) intermolecular potential 
ranges. The obtained results agree well with the ones of computer simulations. 
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1. Introduction 

Binary mixtures in contrast to their constituent components can exhibit three different 
types of two-phase equihbria: vapour-hquid, hquid-hquid and gas-gas The 
possibihty of the reahzation of these phenomena and their priority depend both on the 
external conditions and microscopic parameters of a mixture. The study of tne influence 
of interparticle interactions on the critical properties of a binary mixture is an interesting 
and actual problem. During the last decade this problem has been intensively studied by 
integral equation methods However, this approach, although reproduces different 

phase diagram types by varying the microscopic parameters, gives only a qualitative 
picture of the phenomena under consideration. On the other hand, the critical properties 
of simple fluids and binary mixtures have been studied recently using Monte Carlo (MC) 
simulations 0]-|]l4l. In the vapour-hquid critical temperature was calculated for the 
symmetrical mixture of hard spheres interacting via the square-well potentials. Thus, 
it is interesting to test a theory using such a simple binary fluid model. 

In the present paper we propose a microscopic approach to the study of the vapour- 
liquid critical point of a symmetrical binary mixture. This approach is based on the 
method of collective variables (CV). This method was developed in []T5[ and appears to 
be successful in describing the second order phase transition of the 3D Ising model 



and the vapour- liquid critical point of a one-component fluid [O. On the basis of this 



approach both universal and non-universal quantities were obtained. 

In |]18[ the CV method with a reference system (RS) was generalized for the case of a 
grand canonical ensemble for a multicomponent continuous system. Using this approach 
the phase diagram of the symmetrical mixture was examined within the framework of 
the Gaussian approximation 

In this paper we determine an explicit form of the effective Ginsburg-Landau- Wilson 
(GLW) Hamiltonian of the symmetrical binary mixture in the vicinity of the vapour- 
liquid critical point. Then we integrate the functional of the grand partition function by 



the use of the layer-by-layer integration method proposed in [|16| for the 3D Ising model. 
As a result of this integration one obtains recursion relations for the coefficients of the 
GLW Hamiltonian. The analysis of these relations yields an equation for T^. Here we 
avoid extensive consideration of the results pertaining to the Ising model and call the 



readers' attention to [16| where this problem was studied in detail. The method which 



we describe here yields the same critical exponents as in (see table 1). 

Table 1. Values of the critical exponents and the ratios of the critical amplitudes for 
the 3D Ising model obtained within the framework of the CV method. 
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1.275 
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The paper is organized as follows. We give a functional representation of the grand 
partition function of a two-component continuous system in section 2 and appendix A. 
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In section 3 we construct the basic density measure (GLW Hamiltonian) with respect to 
the CV which include a variable corresponding to the order parameter. In this section 
we also present the basic ideas of the method of the partition function integration in 
the vicinity of the critical point. In section 4 we apply our formalism to calculating the 
critical characteristics (temperature and density) of the binary square-well symmetrical 
mixture. The obtained results are discussed and compared with the MC simulation data 
reported recently by N.B. Wilding jMl and E.de Miquel [|l^ . 



2. Functional representation of the grand partition function of a binary 
mixture 

Let us consider a binary fluid mixture consisting of Na particles of species "a" and Nh 
particles of species "b". The system is in volume V at temperature T. Let us assume 
that an interaction in the system has a pairwise additive character. The interaction 
potential between 7 particle at r, and 6 particle at r^ can be expressed as a sum of two 
terms: 

U^s{\ri - rj\) = ^^s{\ri - r,|) + ^^^dr^ - r,|), (1) 

where "^^slr) is a potential of a short-range repulsion and $7,5 (r) is an attractive part 
of the potential which dominates at large distances. 

A functional of the grand partition function (GPF) of the binary homogeneous 
system in the CV method with a RS can be represented as a product of two factors (see 
Appendix A): 



(2) 



where Sq is the GPF of the RS which we suppose to be known. Si is the part of the 
GPF which is written in the CV space: 



{dp){dc) exp + Co- ^ ^ [^^(fc)pkP-k + 



+ 2i/(fc)pkCk + lV(fc)ckC_k]} J(p,c). (3) 

Chemical potentials /i^ = -^{lA + lA) ^'^d = -^{lA — lA) ^^re determined from the 
conditions: 



dlnSi 
dPfif 
dlriEi 
dPpi 



< Na> + < Nb >=< N > (4a) 

< AT, > - < iV, > . (46) 



Functions V{k),W{k) and U{k) are combinations of Fourier transforms of the initial 
interaction potentials ^^s{k): 
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V{k) 
U{k) 
W{k) 



' 2 ' 
^ 2 ' 



[aaa{k) + abb{k) + 2a„;,(A;)] , 

[«aa(/i;) + Olbb{k) - 2aab{k)] . 



(5) 



J(p, c) = j {dv){duj) exp < i27r ^(cJkPk + t'kCk) + 



n>l i„>0 



(6) 



is a transition Jacobian to the CV pk, Ck averaged on the RS, variables ci^k, t'k are 
conjugated to variables pk, Ck, respectively. 



(-i27r)" 



n/2 



X 



X t'ki • • • t'k.^t^k.^ + i • • • t^k„(5ki+...+k„- 



(7) 



Index in indicates the number of variables i^k in the cumulant expansion (^. Cumulants 
A4n"''(0, 0) are linear combinations of the initial cumulants A^^^...^„(0, 0) (7i = 
a, 6) (see Appendix B). 

We consider a symmetrical binary fluid mixture (SBFM), i.e. a system in which 
the two pure components "a" and "b" are identical and only interactions between the 
particles of dissimilar species differ. Notwithstanding its simplicity, the SBFM exhibits 
all the three types of two-phase equilibria which are observed in real binary fluids, 
namely: vapour-liquid, liquid-liquid and gas-gas equilibria. For the SBFM lJ{k) = in 
(§) and there are only terms with even indices z„ in the cumulant expansion (0) ||20|| . 



3. The method 



As it was already shown [^, the phase diagram of the SBFM consists of three ranges 
(see figure 1): (1) gas-gas separation and vapour-hquid phase transitions; (2) vapour- 
liquid and liquid-liquid phase transitions; (3) vapour-liquid phase transition only. 

The order of priority of the vapour-liquid and separation phase transitions depends 
on both the external conditions and the microscopic properties of the system. There 
exist two critical temperature branches in such a system: branch {T^~^) connected 
with the variable po and branch {T^'^^) connected with the variable cq |2^. All 
the thermodynamic functions of the SBFM are symmetrical with respect to the 
concentration x=0.5 and have an extremum at this point The concentration 

x=0.5 is a critical one for this model mixture. 
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Figure 1. Three phase regions of the symmetrical mixture depending on the 
microscopic parameters: (1) gas-gas and vapour-liquid phase transitions {T-!~^ > 
T"~'); (2) vapour-hquid and liquid-liquid phase transitions (T"~' > Tfr'-); (3) vapour- 
liquid phase transition only. S+ is the structure factor of the reference system. 



In this paper we consider a symmetrical fluid mixture within range 2 or 3 of its 
phase diagram. 

In the case of the SBFM the variables po and Cq are connected with the order 
parameters for the vapour-liquid and separation phase transitions, respectively |20 



This fact allows us to separate CV pk and Ck into essential and non-essential ones 
depending on the phase transition considered. Since we are interested in the vapour- 
liquid critical point, we can consider CV Ck (and z/k) to be non-essential (CV Ck do not 
contain a variable connected with the order parameter, the coefficients standing at the 
second power of Ck (and z/k) are negative) and we can integrate over Ck (and z/k) with 
the Gaussian density measure. In respect to CV pk it is necessary to construct the basic 
density measure taking into account higher powers of pk (we shall consider a p^ model). 

As a result of integrating in (^-(0) over variables Ck (and z/k) we obtain for the 
GPF: 



where 



j {dp)exp |/5p^po - ^Xl^(^)^kP-k| J{p), 



(8) 



n , , p) 



k A/l + /3iy(A;)A^f (0)/2 



J(p) = /"(rf^)exp|z27r J]^kPk + Xl^^T^ (^) ^ 

Ik n>l ^' ^ / 

X Yl ■^niO)uJk, . . .Wk„5ki+...+k„ > , (10) 
ki...k„ J 

Mn{0) = M^^\0) + AMn. (11) 
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AA^„ are the corrections obtained as the result of integration over variables Ck: 

M'i\0) 1 



AMi 



12 < > 



12 <N > 



72 <N>^ ^ 

k 

+ ^^^^^^^^7^ E^(|k|)^(|ki + k|)^(|k2 - k|), 
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AM. 



+ 



96 



<N> 



<N> 



X 



where 



X ^(|k2-k|)^(|k3 + ki + k|), 

(3<N> W{k) 



In figure 2 the typical behaviour of the potenial V^(A;)/|V(0)| is shown. 




(12) 



(13) 



Figure 2. The behaviour of the Fourier transform V^(fc)/|V^(0)| of the attractive part 
of the interaction potential V(r). 

Let US further assume that V{k) = at |k| > B. Then, integration in (P) over pk 
with |k| > i? leads to 5 - functions and the expression for S contains only the sums over 
k with Ikl < B. 



We consider a set of k vectors, |k| < B, as corresponding to the sites of a reciprocal 
lattice conjugated to a certain block lattice {r^} with A''^ block sites in the periodicity 
volume V: 

y y {Baf<N> 
<Nb>=-^ = = ^ '—^ , 14 

1] = ^po"^ is fraction density. Therefore, one may consider quantity B as the size of the 
first Brillouin zone of this block lattice. 
The shift 

= + A5k, 
Pk = + ^i6k, 

where 



Ml = Mi{0) - 



27rA^4(0) 
A<2(0)A^3(0) MliO) 



~ ^(V^)^ ' ^ ~ ■'-'■■■''^)' transforms ^ into a form containing terms A^i(O), 
A^2(0) and Al4(0) only (the primes on pk and cuk are omitted for clarity): 



J exp I p*po - f 5^ V{k)pkp-k- 



k<B 

+ i2n J2 ^kPk - ^^2(0) J2 ^k^-k - 41 . 1-^4(0)1 X 

k<B k<B ■ 



X ^ t^kit^k2t^k3^k4^ki+...+k4 I" ('^^^)^''('^P)^^'', (15) 

ki...k4<B J 



where 



7Wi(0)7W3(0) 7^2(0)7^2(0) MjiO) 



M^iO) 2MliO) SMliO) j ' 

M2iO) = M2iO) ^'^^^ 



2A^4(0) 

7^4(0) = <Nb> M^iO). (16) 

Expression (p!5| ) for S corresponds to the Ising model in the external field (oi — Ppi) 
with the only difference: cumulants A^2(0), Ai4{0) are functions of the fraction density 
7], temperature T and parameters of the attractive interaction ^^^(/c). 
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After integration over cuk we obtain the following form for the GPF: 



1 <Nb> 



X I eME^ipWpf^ . (17) 



Here 



a4 



X ^ Pki • • • Pk45ki+...+k4 + • • • , (18) 

ki...k4<B 

d2{k) =a2 + pV{k), a2 = J J K{x), 

V 1-^4(0)1 

1X4(0)1 



where 




U{l,x)/U{0,x), 
SK^ix) + 2xK{x) - 2, 

7^2(0). (20) 



U{a,x) is a parabolic cylinder function pSf. Expressions ([T7|)-(pOD have the same forms 



as similar expressions for a one-component system obtained in |T^. This coincidence is 
achieved due to the symmetry of the model under consideration. Ei{p) is the Ginzburg- 
Landau- Wilson Hamiltonian for the SBFM in the vicinity of the vapour-liquid critical 
point. 

In order to integrate the GPF (p^-(pOD over p^ and determine the critical 
temperature we use the method developed in [p!6| , p4| for the Ising model. The essence of 



the method consists in subseqent integration over the layers of the CV space, beginning 
from pk which correspond to short-wave fluctuations. Variations of the coefficients of 
E/i{p) as the result of integration over pk in n subsequent layers of CV phase space are 
described by the recursion formulae derived in |TB[. For the case T > in the interval 
[0,5] there exist three characteristic regions [|T^. The first region Bm^ < k < B 
corresponds to the strongly correlated fiuctuations pk, their density measure is non- 
Gaussian. The procedure based on the renormalization group symmetry is valid here. 
This is the region of the critical regime (CR). The second region < A; < B^^ is related 
to the fiuctuation distributed according to the Gaussian density measure. This is the 
limiting Gaussian regime (LGR). 
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The third region consists of the point k = 0. The variable po is a macroscopic one 
and corresponds to the fluctuations of the particle density in the "external field" fi*. 

We integrate (|1^) according to the following scheme [0. The region (0, -B) is 
divided into the intervals (Si, 5), (Sj+i, -Bj), where i?„ = B / (S is a division 
parameter). Each interval corresponds to a layer of subscripts k in the Brillouin zone and 
each layer of subscripts k - to a layer in the phase space pk- Integrating gradually over 
the layers we get a block lattice sequence with an appropriately growing block period and 
with the Hamiltonian corresponding to each block. Each Hamiltonian is characterized by 
the coefficients d2, a^, ^2^'', a^^^; d^2 \ \ stc. For the sequence of the block Hamiltonians 
I 5 '^4"'' I the renormalization group symmetry holds and the fixed point is of a saddle 
type. Because the explicit expressions for the initial values of coefficients d{k) and 04 
are given (see (p!^)-([20|)), the solutions of the renormalization group type are functions 
of microscopic parameters, density and temperature. 

The CR takes place for all the variables pk at the critical point. Therefore, the 
critical temperature can be determined from the solution of recurrent equations. Using 
the result obtained in we derive the formula 

A(/?V(0))2 + 5(/?r(0))+D = 0, (21) 

whre 

B = -02, 

D = a4i?(°Vv^- 

/o, (fo are coordinates of the reduced fixed point, R^^^ is a universal function of parameter 
S. The optimal value of S is 3.4252 and the values of fo,^po, R^^^ corresponding to it are 
taken from From the condition p* = we obtain the second equation ||17|| : 



MsiO) = 0, (22) 
which allows us to determine the critical density of the system. 

4. Results and discussions 

In this section we present our results for the vapour-liquid critical point of symmetrical 
mixtures, using the method proposed above. These results are compared with those 
previously obtained by Monte-Carlo simulations [1^, [1^. 

The system under study is a symmetrical hard sphere square-well binary mixture. 
The interaction potential between the particles is given by: 

00, if r < a 

U-ysir) = ^ -e^5, if a < r < Xa , 
0, if r> Xa 

where cr is a hard sphere diameter, A is a range of the potential, and well-depth 
of the interaction between the particles of types 7 and 6. For a symmetrical mixture 
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eaa = ^bb = ^ 1^ ^ab- In OUT formalism a completely analytical treatment for general A is 
possible. But in this paper we choose A = 1.5 and A = 2.0, because for these values of 
A the MC results are available. 

We split the potential U^s{r) into short- and longe-range parts using the Weeks- 
Chandler- Andersen partition [^. As a result, we have: 

%sir) = I : - ! . (23) 




*•,<('') = •! 7" zi t: ■ (24) 

In our case the RS is a hard sphere system with the diameter a (see (PBD). The 



Fourier transform of function (24) has the form: 



~ 3 

^-ysik) = $^g(0) ^^ Ax cos(Ax) + sm(Ax)], 

where 

X = ka, 

Cumulants Mt\o,:;0) are calculated according to the formulae given in 
Appendix B. For 6*2(0) the Percus-Yevick approximation is used: 

^ ' (l + 2r/)2 

The solutions of equations (PT|)-(P^ are found numerically using a self-consistent 
procedure by means of which the dependences of the coefficients 02 and 04 (as well as 
cumulants A^„(0)) on (3^ are taken into consideration. 

The vapour-liquid critical temperatures Tc (T^ = ksT/e) versus the microscopic 
parameter r (r = Cab/e is a dissimilar interaction strength) are shown for A = 1.5 and 
A = 2.0 in figure 3. 

In figure 4 we demonstrate the dependence of the critical density rj^ on r for A=1.5 
and A=2.0. We compare our results with those obtained from MC simulations: for 
A = 1.5 and r=0.72 we have = 1.055, while the MC simulations give = 1.06(1) 



1^ , for A = 2.0 and r=1.0 (the case of r=1.0 corresponds to a one-component system) 
we obtain Tc = 2.753 and r/c = 0.129, while the simulations give Tc = 2.684(51) and 
r]c = 0.123(43) ^ 



5. Conclusions 



In this paper we propose a method for the study of the vapour-liquid critical point 
of a symmetrical binary mixture depending on its microscopic properties. Using this 
method we calculated the critical temperature and critical density of the symmetrical 
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Figure 3. The vapour-liquid critical temperature as a function of the microscopic 
parameter r at A = 1.5 (left) and A = 2.0 (right). 




Figure 4. The critical density as a function of the microscopic parameter r. 



hard sphere square-well fluid versus the microscopic parameter r which is the measure 
of the dissimilar interaction in the system. Our results agree with those obtained by 
MC simulations. We can improve our results in the following ways: 1) taking into 
consideration the region of k with |k| > i? (see figure 2); 2) using a higher approximation 
than the one. 

Appendix A 



A grand partition function of a two-component fluid system in the CV representation 
with a RS can be written as in HTSl: 
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where 



oo oo 



Na=0 Nt,=0 7=a 



NJ 



(dT) exp 



7,(5=a,6 i,j 



is a grand partition function of the RS; (3 — /cs- is the Boltzman constant, T is 
temperature; (dF) = Ha 6 ^^^7' ^^^^^ = drldf^ . . . drjf is a volume element of the 
configurational space of the 7-th species; /Iq is a chemical potential of the 7-th species 
in the RS. 

The part of the grand partition function which is defined in the CV phase space 
has the form of a functional integral: 

Si= /(c?p)exp[/3^/i>o,7-^ Yl Xl"^^(^)^k,7P-M]^(Pa,Pfe)-(A.l) 



7,5=0,6 k 



Here, 



1) iij is a part of the chemical potential of the 7-th species 

= /^7 - /^O + ^ X] «77(^)' 

and is determined from the equation 
9 In Si 



/i^ is a full chemical potential of the 7-th species; a^s{k) = ■^^^^(A;) ; < A^^ > is an 
average number of the 7-th species particles. 

2) pk,7 = Pk 7 ~ ^Pk 7(7 ~ ^) collective variables of the 7-th species, the indices 
c and s denote the real part and the coefficient at the imaginary part of pk,7; Pk,7 

^ describe the value of the k-th fluctuation mode of the number of the 7-th species 
particles. Each pj^^ and takes all the real values from —00 to +00. {dp) is a volume 
element of the CV phase space: 

(dp) ^Yldpo^Yl dpl^^dpi,^. 

7 k^^O 

The prime means that the product over k is performed in the upper semi-space; 

3) J{pa,Pb) is a transition Jacobian to the CV averaged on the RS: 



b 

J{Pa,Pb) = / (rfz/) JJexp 

7=a 



i27T ^ t'k,7Pk,- 



exp 



L 



.n>l 



ni 



^71 ■••7n (ki, ■ ■ ■ , k„)i/ki,T-i . . . i/k„,^„ 

7i...7„ ki...k„ 

where variables z/k,^ are conjugated to CV pk,^- Af^i...^^ (ki, .... k„) is the n-th cumulant 
connected with >S'^i...7„(/i;i, • • • , kn), the n-particlc partial structure factor of the RS, by 
means of the relation 

^7l-7n(kl, ■ ■ ■ ,kn) = \/^7i • • • ^7n'S'7i...7„(fcl, • • • , /i;n)^ki+-+k„ , 
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where 5kiH hk„ is a Kronecker symbol. 

In general, the dependence of Myj...^^(ki, . . . , k„) on wave vectors 
ki, . . . , kn is complicated. Hereafter we shall replace M^^...^^(ki, . . . , k„) by their values 
in the long- wave length hmit M^j .,^^(0, . . . , 0); 

4) (j)^s{k) is a Fourier transform of the attractive potential (j)^s{f)- Function 
(p-fsik) satisfies the following requirements: (p^s{k) is negative for small values of k 
and Zzmk_,oo075(A;) = 0. 

We pass in ( |A.1| ) to CV pk and Ck (according to and 7k) by means of the 
orthogonal linear transformation: 

Pk = ^(Pk,a + Pk,fe), Ck = — (Pk,a - Pk,fe), 

UJk = -^{^Ka + t'k.fc), t'k = -^{^Ka - t'k.b)- (A. 2) 

Now Pk and Ck are connected with the total density fluctuation modes and the relative 
density (or concentration) fluctuation modes, respectively. 
As a result, for Si we obtain formulae (^-(0). 



Appendix B 



Cumulants Ain"\o) with n < A are expressed in terms of the initial cumulants 
'^7i...7n(0, •••,0) (71, ...,7„ = a, 6) as follows 0]: 

0) = MaiO) + MbiO) =<N> 
0) = ^„(0) - MbiO) =<Na>-<N,> 

0)=MaaiO)+MbbiO)+2MatiO) 
0) = MaaiO) - MbbiO) 

0) = MaaaiO) + MbbbiO) + 3[MaabiO) + Mabb{0)] 
0) = MaaaiO) - MbbbiO) + MaabiO) " MabbiO) 
0) = MaaaiO) + MbbbiO) " MaabiO) " MabbiO) 
0) = MaaaiO) " MbbbiO) - 3[MaabiO) " MabbiO)] 
0) = MaaaaiO) + MbbbbiO) + 

+ ^MaaabiO) + MabbbiO)] + GMaabbiO) 
0) = MaaaaiO) " MbbbbiO) + 2[MaaabiO) - MabbbiO)] 
0) = MaaaaiO) + MbbbbiO) - 2MaabbiO) 
0) = MaaaaiO) " MbbbbiO) " 2[MaaabiO) " MabbbiO)] 
0) = MaaaaiO) + MbbbbiO) - 

- ^[MaaabiO) + MabbbiO)] + QMaabbiO) (B.l) 

The same expressions hold at kj 7^ 0. 



(o; 



M'l 
M^^ 
Mf 

mS 

Mf 
Mf 

(2 



M 

Mf 

(o; 



M 



(1 



M 

Mf 
M^' 
M^' 
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The nth cumulant A4n"\o) with i„ = is connected with the nth structure factor 
of the one- component system S'„(0) [pO[] : 

M^^\0) =<N>Sn{0). 

Structure factors 5'„(0)(n > 2) can be obtained from 82(0) by means of a chain of 
equations for correlation functions Cumulants with z„ 7^ can be exspressed in 

terms of A4.n\o) (see formulae (4.8) in |^0|| ) . 
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Figure captions 



Figure 1. Three phase regions of the symmetrical mixture depending on the 
microscopic parameters: (1) gas-gas and vapour-liquid phase transitions (TJ'"^ > 
T^~^); (2) vapour-hquid and liquid-liquid phase transitions (TJ~' > T^~'); (3) vapour- 
liquid phase transition only. S+ is the structure factor of the reference system. 



Figure 2. The behaviour of the Fourier transform y(fc)/|y(0)| of the attractive part 
of the interaction potential V(r). 



Figure 3. The vapour-liquid critical temperature as a function of the microscopic 
parameter r at A = 1.5 (left) and A = 2.0 (right). 



Figure 4. The critical density as a function of the microscopic parameter r. 



